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We confront for the first time the widely-held belief that combined event-by-event information 
from quark gluon plasma signals can reduce the ambiguity of the individual signals. We illustrate 
specifically how the measured antiproton yield combined with the information from pion-pion HBT 
correlations can be used to identify novel event classes. 
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I. INTRODUCTION 

Signals of novel phenomena from hadron production in 
relativistic heavy ion experiments have proven ambiguous 
at the Brookhaven AGS and the CERN SPS, primarily 
because descriptions of the 'ordinary' hadron production 
mechanisms are under-constrained. It has long been con- 
jectured - but not yet demonstrated - that the added 
information from the correlation of distinct signals can 
reduce the ambiguity of the individual signals. In partic- 
ular, the measurement of such correlations is one of the 
driving principles behind the PHENIX and STAR exper- 
iments at RHIC JD. However, to describe the correlated 
information that these experiments yield, phenomenolo- 
gists must introduce new unconstrained parameters. 

We propose that added information can in fact be 
gained by studying antiproton production in conjunction 
with pion interferometry, HBT, on an event-by-event ba- 
sis. Various authors have suggested that both antiproton 
production || and HBT radius parameters J3| can change 
abruptly and dramatically if quark gluon plasma forms in 
a heavy ion collision. However, theoretical uncertainties 
|Q-g| and experimental difficulties M have made results 
difficult to interpret in that context. In this paper, we 
demonstrate the utility of covariance measurements by 
constructing a plausible scenario in which a measurement 
of the covariance of these observables can resolve ambi- 
guity in individual HBT and antiproton measurements. 
After explaining the basic scenario, we discuss how an- 
tiproton production and HBT radii become correlated in 
the framework of a thermal model of hadron production. 
We then develop a Monte Carlo code to provide a realistic 
simulation of our scenario in the context of STAR. 

Measurements of the covariance of distinct signals are 
useful when event averaging hides otherwise-strong sig- 
natures of a new event class. Suppose that there are two 
event classes - "plasma" and "hadronic." Further assume 
that the mean antiproton rapidity density N = dn/dy is 
different in each class, so that the signal is truly strong. 
Averaging over events yields 



plasma and hadronic classes respectively. The event- 
averaged (N) is a smooth function of centrality and beam 
energy, because / is a continuous function of the collision 
geometry and energy deposition. Smooth data sets, such 
as those in figs, la and lb, are subject to broad inter- 
pretation. 
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FIG. 1. Antiproton multiplicity (a), pion HBT transverse 
radius (b) and their covariance (c) as functions of the charged 
particle multiplicity. Results are obtained from simulations 
below for the scenario I, cf. eq. (Eq). Multiplicities are com- 
puted for the STAR acceptance. 

The correlation of N with the HBT transverse radius 
Rt can betray the existence of the new event class. This 
correlation is characterized by the covariance: 



(N) =J2N t = fN(q) + (1 - f)N(h), 



(1) 



C nr = J2(Ni - (N)){R T i - (Rt)), 



(2) 



where the / is the probability that the i th collision forms 
a plasma and N(g, h) are the average values for the 



where Rt is the transverse radius measured in event-by- 
event pion interferometry. We find 



AC„, = /(I - f){N(q) - N(h)}{R T (q) - R T (h)}, (3) 

where C = fC(q) + (1 - f)C(h) + AC and the covari- 
ance for each class is C(q,h). The term AC depends on 
the hidden shift in N and Rt- In the best of all possi- 
ble worlds, / will change from zero to unity for a single 
target-projectile combination as the impact parameter b 
is varied. This variation can introduce a peak as shown in 
fig. lc, provided that both class-averages Rr(q) and N(q) 
truly exceed the hadronic values, as predicted. More 
likely, one may come upon a low impact-parameter re- 
gion where / begins to rise, or a peripheral region where 
/ begins to fall. 

We point out that the variances cr-^L and <j\ can 
exhibit similar behavior in the presence of two event 

2 

classes. For antibaryons, we find (AN ) = fazjq) + 
(1 - f)a^(h) + A4, where 



Aa w 2 = /(l - f){N(q) - N(h)Y 



(4) 



q.h 



and a— are the standard deviations for each class; the 
result for Rj- is similar. The variance a^- can be par- 
ticularly interesting in the context of recent work by 
Stephanov, Rajagopal and Shuryak ||; we address that 
point elsewhere PJlO[]. 

A measurement of (|2|) is built upon an event-by-event 
analysis on which we now comment. We expect any- 
where from 50 to 100 antiprotons per event in the STAR 
acceptance, and possibly more pdj . Such a number is 
perfectly adequate for event-by-event analysis. Since our 
intent is to use the antiproton yield as a proxy for an an- 
tibaryon measurement, we are not concerned if the mea- 
surements contain a contribution from antilambdas, as 
found at the AGS @. On the other hand, the HBT part 
of the measurement is challenging because it involves an 
event-by-event two-pion correlation analysis. In an HBT 
analysis, one compares the measured identical-pion cor- 
relation function C^n to a Gaussian parameterization. 
Radius parameters are typically obtained by a three di- 
mensional fit to high-statistics data (see B for details). 
The hadronization of the plasma affects these spacetime- 
dependent variables directly by delaying the freezeout of 
the system (p). While such a three-dimensional analysis 
may not be practical for single RHIC events, the feasibil- 
ity of one dimensional analyses has been studied by the 
STAR experiment fll2|l . One can extract the transverse 
radius by comparing 



CV.(<ZT) = l + e-^ 



(5) 



to data, where qr is the difference in the pions' trans- 
verse momenta. Alternatively, Heinz and Wiedemann 
suggest that it may be easier to extract event-by-event 
parameters using R^ 2 oc J d 2 qTq^{C(qT) — 1) (l3). The 
difference between these definitions is immaterial to our 
discussion. However, we point out that Rt, which essen- 
tially measures the transverse system size, is not an ideal 



plasma probe. Nevertheless, we find a strong change in 
C nr for a 50% increase in Rr(q) due to the transverse 
growth of the longer-lived plasma system. 



II. FLUCTUATIONS NEAR EQUILIBRIUM 

To establish the mechanisms that drive the correlations 
and fluctuations in antibaryon production and pion HBT, 
we employ the idealized but standard Bjorken hydrody- 
namic framework that describes particle production near 
rapidity y = 0. We assume that matter in this region is 
in local thermal equilibrium, with an average local tem- 
perature T, entropy density s and net baryon density ps 
that vary only with proper time r. The transverse area 
A is initially determined by the overlap of the colliding 
nuclei. One can define a comoving volume for matter in 
the central region, V = S/s = At. This volume grows 
from a formation time To to freezeout at tf, so that the 
entropy per unit rapidity for all hadrons S is r indepen- 
dent. The total rapidity density of hadrons AT tot is nearly 
constant, because A^ot oc S for a system dominated by 
light hadron species with masses <C T. Baryon current 
conservation implies that the net baryon rapidity den- 
sity Nb = dns/dy = Apb(t)t is independent of proper 
time. Mean rapidity densities of individual species, such 
as antibaryons N and baryons N, vary with r. Observe 
that this model is distinct from global equilibrium mod- 
els used by many groups. The key distinction is that 
here, we focus on thermodynamic quantities local to the 
central region. 

In local equilibrium, the number of antibaryons at 
midrapidity and the source size fluctuate depending on 
the fluctuation of the state variables in this region. We 
consider an ensemble of collisions at fixed impact parame- 
ter in which V, T and Nb fluctuate |I3. We assume that 
these variables are independently established near y = 
in the preequilibrium evolution, so that they are statisti- 
cally independent. The comoving volume and tempera- 
ture fluctuate because the initial number of mesons and 
the initial energy per meson vary from event to event. 
The net baryon number fluctuates because the central 
region can exchange baryon current with the rest of the 
collision volume for fixed V and T, i.e. the central re- 
gion is held at constant baryon chemical potential. We 
remark that this ensemble differs marginally from the 
familiar Grand Canonical ensemble, in which the total 
hadron energy fluctuates and T is held fixed. The dif- 
ferences are negligible in the regime Nb «C A to t that we 
consider. 

The fluctuations of a quantity X are characterized by 
the variance a\ = (AX 2 ) for AX = X - (X). The 
fluctuations of V satisfy 
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(6) 



where A to t is the total number of hadrons. Thermal fluc- 
tuations satisfy p4[: 
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for a system of a heat capacity C v , where the second 
equality follows for an ideal gas of (mostly) massless 
hadrons in which N B <C N to t and Fermi and Bose statis- 
tics can be neglected. The variance of the net baryon 
number at constant temperature and volume is 



(AN 



TdN B /d^ B 



(8) 



by straightforward extension of rcf. pi . For an ideal gas, 
N = Vf(T) exp(- MB /T) and N B = 2Vf(T) sinh( MB /T) 
(we need not specify /), implying that 



a B = N + N, 



(9) 



neglecting small corrections from Fermi statistics. 

We comment that Stephanov, Rajagopal and Shuryak 
pf have discussed thermal fluctuations of observables at 
a critical end point, focusing on the effect of the diver- 
gence of Cy on (RV We point out that critical fluctu- 
ations also affect (g), because 8N B /dfi B = —d 2 Q/dfi 2 
diverges at that point, where fl is the free energy. The 
antibaryon fluctuations that we discuss shortly would re- 
flect this critical behavior, c.f. ( |l2| ) below. Here, we as- 
sume that T and [i B are sufficiently far from that point 
that these effects are negligible. 

The extent to which the number of antibaryons fluc- 
tuates depends on whether or not chemical equilibrium 
holds with respect to the reactions that change the num- 
ber of baryon-antibaryon pairs, e.g. NN ^ mesons. We 
do not expect chemical equilibrium for these processes 
unless a plasma is formed |lq |. Suppose then that the 
numbers of baryons and antibaryons are individually con- 
served. For constant T and V, ref. JlJ] implies that the 
antibaryons essentially follow Poisson statistics, 



(AN 



v =N, 



no chemical equilibrium (10) 



with a similar expression for protons. Observe that HI- 
JING and similar event generators exhibit exactly the 
same Poisson-like behavior (such behavior is built in to 
these models). 

Chemical equilibrium couples fluctuations of the 
baryons to those of the antibaryons, reducing the relative 
variance compared to (pOJ) . In this case, a small change 
in the number of antibaryons requires the net baryon 
number change: 



AN=(^?^ \AN B . 



\dN B /dfi B 

The variance in antibaryons is then 



(11) 



AV /; -T^j^ = N(N + N)-\ (12) 



where we have used (|9j). Note that we recover (|10|) in 
the limit N <C N, since baryon conservation keeps N 
constant. 



Volume and thermal fluctuations can also cause the 
number of antibaryons to fluctuate. We write: 
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(13) 



The third contribution represents the constant T and V 
results discussed earlier. The second term, 



dN 
~dT 



dN 



/ dN \ dN B /dT 



\d^ B ) ON B lm 



(14) 



allows thermal fluctuations to change the antibaryon 
population by making pairs in chemical equilibrium; it 
strictly vanishes otherwise. We will see that the covari- 
ance (||) arises from the first term. 

For an ideal gas in chemical equilibrium, 



AN 



AV 
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2eN AT 
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N+N T N+N 



(15) 



where e = Ej^/NT w m/T + 3/2 is the energy per an- 
tibaryon per unit temperature. We then obtain: 
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W (I 6 ) 



The third term, i.e. (|l2|), dominates the antibaryon fluc- 
tuations when JVtot ^> N + N, as eqs. (gM suggest. For 
N 3> N we obtain the generalization of (Q) , 

c— a 2 1 

— % = -j^ + =. no chemical equilibrium (17) 
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The underlying volume fluctuations follow (|6|) provided 
that thermal equilibrium holds. 

We expect volume, thermal and baryon-number contri- 
butions to be generally comparable at RHIC in thermal 
equilibrium. The sum of the volume and thermal terms 
is roughly 1 + e 2 /12 sw 7 for T = 140 MeV and m = 
938 MeV, implying a contribution to (a w /N) 2 of - 0.7% 
for TVtot ~ 10 3 as expected in Au+Au collisions. Baryon 
density fluctuations contribute ~ 1.3% to the variance 
for N w N w 40, yielding a total (<r w /N) 2 - 2%. Ob- 
serve that either the chemical nonequilibrium result (nfl) 
or HIJING would give a similar value ~ 2.5%. 

We expect the covariance (0) to be determined primar- 
ily by volume fluctuations, as we see by computing the 
related covariance of N and V, 



r = 



(ANAV) = Nay /V. 



(18) 



The large contributions to the antibaryon variance (|lq ) 
from net baryon number and temperature fluctuations 
are absent in (18). A very important consequence of 
this result is that this expression is valid whether or not 
chemical equilibrium holds. 



To relate the volume fluctuations to fluctuations in the 
HBT radius, we observe that variations in HBT trans- 
verse radius are driven primarily by the geometrical ef- 
fective radius ~ A 1 / 2 and, secondarily, by flow [pi. By 
analogy with (fl3), we make the phenomenological ansatz: 



AR T /R T = kAV/2V + AAT/T, 



(19) 



where the parameters k — 2d\ogRT/dlogV and A = 
d\og Rt I d\ogT generally must be determined by hydro- 
dynamic calculations. We take the relative concentra- 
tions N/Ntot and N/N to t to be small enough that the 
baryons have a negligible effect on the flow at RHIC. 
The covariance of Rt and N is then 



C nr = (ANARt) = NR T k 



"2V 2 



+ e\ 



T 2 



while the fluctuations in Rt itself satisfy 
a\ T jR T = n 2 a 2 v /W 2 + A 2 4/T 2 



(20) 



(21) 



In contrast to (jlq), there is no first-principles justification 
for (g(j, ^l|) - flow effects are outside the reach of our 
thermodynamic approach. 

In the absence of transverse flow, Rt is essentially the 
geometric transverse radius g, so that flow contribu- 
tions to ( p0| ) and (EfT) must vanish. In keeping with our 
Bjorken-like scenario, we therefore take 



1 and A = 0. 



(22) 



To see that this approximation is reasonable, we use a pa- 
rameterization of the Yano-Koonin-Podgoretskii radius 
Rt(Mt), eq. (47) of ref. H, which is closely related to our 
Rt- For that parameterization, n is strictly unity. We 
estimate A ~ 0.3 for a plausible mean transverse fluid ra- 
pidity of r)f = 0.6. This result is in keeping with findings 
at SPS energy, where NA49 finds that Rt(Mt) decreases 
by ~ 20% as wit is varied from zero to sw 0.45 MeV, cor- 
responding to A ~ 0.2. One can improve our estimate by 
introducing stochastic noise and dissipation into a hydro- 
dynamic model and studying the consequent fluctuation 
of R T - 

We comment that by using the comoving volume V 
for the Bjorken scaling expansion, we can describe the 
local equilibrium fluctuations in our evolving system as 
if the system were stationary. However, for more general 
local equilibrium systems flowing in three dimensions, 
i.e., real heavy ion systems, we must obtain the vari- 
ance and covariance from the local hydrodynamic fluctu- 
ations 16 1 . To see how inhomogeneity can affect our esti- 
mates, consider temperature fluctuations that locally sat- 
isfy (AT(x)AT{x')) = [T{x) 2 /n(x)c v (x)]5(x-x'), where 
n is the total hadron density and c v the specific heat. 
Let us assume that longitudinal Bjorken scaling holds 
but that the system is inhomogeneous in the transverse 
direction r - the simplest possible extension of our ideal- 
ized model. Experimental measurements of the average 



transverse momentum roughly yield the volume average 
T ex J d 2 rn(r)T(r) in each event |J . We can then use this 
volume average to construct the event average (AT 2 ). 
Using the above correlation function, we see that this 
average satisfies (|7|), since local equilibrium implies that 
T(r) 2 /c v (r) scales as n _1 , neglecting corrections from 
Bose and Fermi statistics. In this case and in the case of 
rapidity densities, we do not expect transverse inhomo- 
geneity to introduce large corrections. However, experi- 
ments do not strictly yield the volume average of T but, 
rather, momentum distributions from which we extract 
an average pt slopes. To treat fluctuations of pt slopes, 
Rt and other such quantities with precision, we must use 
hydrodynamic or transport models that include noise and 
dissipation that respect the fluctuation-dissipation theo- 
rem [Q6l. No such models exist. 



III. FLUCTUATIONS IN ION COLLISIONS 

To describe the two-event-class scenario outlined in the 
introduction, we now apply these general results to de- 
velop a Monte Carlo event generator for correlated sig- 
nals in RHIC collisions. With this generator we can ac- 
count for additional sources of fluctuations in heavy ion 
collisions as described below. Schematically, we gener- 
ate events as follows. For each event class, the average 
N, Rt and total rapidity density iVtot are determined 
at each impact parameter by the collision geometry us- 
ing the prescription that we describe shortly. We then 
generate values of Ni,RTi and the charged-particle mul- 
tiplicity for each event in accord with the variances and 
covariance calculated earlier. As possible in the STAR 
experiment, we use the charged particle multiplicity to 
determine the centrality of each event. We then compute 
the covariance as a function of multiplicity for two illus- 
trative two-class scenarios. In addition, we present the 
distribution of simulated events to eliminate some of the 
abstraction that often obfuscates statistical analyses. 

In discussing thermodynamic fluctuations, we have so 
far assumed that collisions occur at a fixed impact pa- 
rameter b, with fluctuations occurring about well defined 
mean values fixed by V, T and Nb- Additional fluc- 
tuations arise in heavy ion experiments, e.g., from the 
imperfect experimental knowledge of the impact param- 
eter. In this work, we incorporate such fluctuations into 
our description using a Monte Carlo framework. How- 
ever, one can understand how these fluctuations come 
about as follows. Events with impact parameters in the 
range from b — Ab/2 to b + Ab/2 yield average numbers 
of antibaryons that differ by « (dN/db)Ab. An experi- 
ment that does not resolve these impact parameters will 



measure a covariance: 






dN\ (dR T 



db 



db 



'b-> 



(23) 



where <7f, is the standard deviation for events in the unre- 
solved range and C%* is the average equilibrium contribu- 



tion for collisions at b. This centrality contribution must 
vanish for central collisions by symmetry. More generally, 
its magnitude depends on how the mean values N and 
Rt vary with collision geometry and, additionally, how 
centrality is measured. For Au+Au collisions in STAR, 
we find that impact parameter fluctuations are typically 
more important than volume fluctuations, but less impor- 
tant than the thermal and baryon-number contributions. 
[At an impact parameter b = 6 fm, we use eqs. (M B9 ) 
to estimate the centrality contribution to the relative co- 
variance Cut/NRt, to be ~ 0.3% compared a thermal 
contribution of ~ 0.05%.] 

To estimate antibaryon production for the purely 
hadronic event class, we observe that kinetic theory 
estimates suggest that it is unlikely that baryons will 
reach chemical equilibrium ]13,nffl. At RHIC, events 
too peripheral to produce a plasma will yield fewer an- 
tibaryons than required by detailed balance for the re- 
actions NN ^ mesons. Chemical equilibration then re- 
quires that meson collisions increase the antibaryon pop- 
ulation, but time scales for those processes greatly exceed 
the relevant dynamic time scales. We can therefore ex- 
pect the hadron fluid to be far from chemical equilibrium, 
although thermal equilibrium may hold. In this case, the 
mean rapidity density of antibaryons is essentially the 
same as its initial value, since antibaryon number is now 
effectively conserved. Moreover, we use (17) and ( jig ) to 
compute the standard deviation crj^(h) and covariance 
C nr (h). Observe that the final state is practically in- 
distinguishable from the entirely-nonequilibrium initial 
state. For Au+Au at RHIC we take: 



N(h) = nhAf(b)/Af(0); 



(24) 



where rih = 40, a value consistent with the range of event 
generator predictions: ffiJING, HIJING/Si? (its succes- 
sor) and RQMD report rapidity densities of 60, 20 and 20 
respectively. We assume the rapidity density varies with 
impact parameter in proportion to the number of partic- 
ipants, M(b), in agreement with ITIJING calculations at 
the 1-2 percent level. 

Antiproton production can be very different in a 
plasma. A variety of mechanisms from chiral restora- 
tion to disoriented chiral condensate formation B can 
enhance production of baryon-antibaryon pairs and fa- 
cilitate equilibration. Model calculations typically yield 
values of the antiproton rapidity densities well in excess 
of event generator estimates. For example, ref. [|ll| pre- 
dicts values of N of about 86, almost three times the 
HIJING level. For an enhancement at that level to be 
hidden in the mean value (pi), / would have to be very 
small. We assume a more conservative 30% enhancement 
over the hadron gas value in a central collision, 



26, 



(25) 



and take the same centrality dependence. Importantly, 
since chemical equilibration is likely in a plasma, we now 



use ( |16[ ) to calculate ojj(q). On the other hand the co- 
variance for chemical equilibrium is still given by ( pL8| ) , so 
that Cnr(q) « C nr (h). 

For our HBT estimate, we take a mean transverse ra- 
dius in the hadronic event to be roughly the geometric 
radius 



R T (h)=r h [A(b)/A(0)} 1/2 , 



(26) 



where r/, — 6 fm and A(b) is the geometrical overlap 
area of the colliding nuclei. For comparison, Hardtke 
and Voloshin J2CJ has used RQMD to find a side radius 
R s sw 5 fm in Au+Au at RHIC and we expect Rt ~ R s - 
We then use (20) and (22) to estimate the fluctuations. 
Pratt and Bertsch have argued that plasma for- 
mation can increase the pion HBT radius parameters, 
e.g., if a nearly- first-order phase transformation leads to 
the dramatic increase of the collision-system lifetime and 
size. However, as noted earlier, this effect is only dimly 
reflected in Rt- We assume a 50% increase over the 
hadronic value, 



9 fm. 



(27) 



We expect the intrinsic uncertainty from (Ell) to be much 
smaller than the experimental uncertainty. 

We assume that STAR will measure the multiplic- 
ity of charged pions to select the centrality of each 
event. We compute the multiplicity for each event as- 
suming a Gaussian distribution with an average value 
N to t(b) = N tot (0)Af(b)/N(0) and a standard deviation 
0tot = y/Ntot consistent with thermal equilibrium. The 
scale iVtot(O) = 2100 is determined by the initial produc- 
tion regardless of event class, in accord with entropy and 
energy conservation; the particular value is taken from 
a HIJING simulation in the STAR acceptance. We then 
generate the antiproton yield and Rt using (6), (7), (12) 
and (18). To each Np *=* N/2 and Rt value we add a 
simulated "experimental" fluctuation. The experimental 
fluctuation in Rt is distributed with AR/R ~ 10%, in 
accord with 12 1, while that of Np assumes an ad hoc 95% 
detection efficiency. 

Importantly, these experimental fluctuations do not af- 
fect the covariance, although they do affect the scatter 
of events. Only correlated uncertainties in the N and 
Rt measurements can alter the covariance. Since these 
two measurements are very different, we do not antic- 
ipate any correlated uncertainty (although real experi- 
ments are very complicated!), as long as the measured N 
and Rt truly come from the same event. 

To illustrate the effect of two event classes, we present 
two ad hoc scenarios for the onset of plasma events in 
RHIC Au+Au collisions. We start with a Pangloss sce- 
nario in which plasma forms with certainty in all but the 
most peripheral collisions. The large fluctuations in stop- 
ping and energy deposition inherent in peripheral colli- 
sions might plausibly result in an event class in which 
plasma does not form. We take the probability f(b) for 
plasma events to be 



/i(6)={l + cxp[(6-6 )/A&]}- 



(28) 



a form that would apply in the best of all possible worlds. 
The plasma fraction increases from zero to one within a 
range ~ 2A6 ~ 1 fm of bo = 6 fm. In figs, la and 
lb we show the average N and Rt derived from 10,000 
simulated events. We see that the average values vary 
smoothly between hadronic and plasma model expecta- 
tions. In fig. lc we present the covariance of Rt from 
10 6 events (STAR will obtain that many events in twelve 
days). Our simulations clearly establish the behavior of 
(pi) within the statistical uncertainties. Observe that the 
tails of the distribution are described mainly by (23), 
with volume fluctuations contributing only at the 10-30% 
level. Further note that in the limit as Ab — > 0, the width 
of the bump in fig. lc tends to the value at t(&o) set by 
the multiplicity distribution. 
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FIG. 2. Antiproton multiplicity (a), pion HBT transverse 
radius (b) and their covariance (c) for scenario II, cf. eq.y29|). 

We now consider a more conservative - but no less ad 
hoc - scenario in which plasma events appear in only a 
fraction of the most central collisions. We take: 



fu(b) = 0.25 [1 - (6/6 ) 2 



(29) 



for b < bo = 3 fm, with / = otherwise. We imagine that 
a region of T > Tc forms in central collisions, and that 
this region grows in size as collisions become more cen- 
tral. The event-averaged Monte Carlo results in fig. 2a 
and 2b show no appreciable difference from hadronic ex- 
pectations, while the covariance in fig. 2c shows a striking 
enhancement. 
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N ch 
FIG. 3. Distribution of events: transverse radius vs. 
charged particle multiplicity for 10,000 simulated events for 
scenario II. 




FIG. 4. Distribution of events: antiproton multiplicity vs. 
charged particle multiplicity for 10,000 simulated events for 
scenario II. 

To appreciate how the event distribution gives rise to 
this behavior, we show the distribution of events in scat- 
ter plots in figs. 3, 4 and 5. While the antiproton results 
in fig. 4 are unconvincing, the HBT events in fig. 3 show 
a hint of a new population in central collisions. The cor- 
related signal from 10,000 events is shown in fig. 5. Here, 
two event classes are clearly distinguishable. The scat- 
ter plot indicates that we can identify the events with 
the largest antiproton yield as anomalous. Armed with 
such information, experimenters can then search for an 
optimal beam energy and target-projectile combination 
to sweep out the coexistence region as in fig. 1. More- 
over, one can introduce an antiproton trigger to other 
experiments to gather a statistically significant sample 
of plasma events, e.g., to study hard probes or three di- 
mensional HBT. We point out that the covariance can 
be measured to arbitrary precision by collecting events 
and, consequently, is more sensitive to the appearance 
of a new event class than is the qualitative scatter plot. 



The result in fig. 2c is computed from f 6 events. 




To demonstrate how correlated signals can help es- 
tablish a new event class, we have taken literally many 
of the predictions of quark gluon plasma formation. 
While other assumptions, e.g., chemical equilibrium, are 
not necessary to our conclusions, the existence of event 
classes with distinctive signals is crucial. The observ- 
ables themselves need not shift abruptly as in a contin- 
uum phase transition, but their values must be markedly 
different in each class - correlating non-signals will not 
help. 
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FIG. 5. Distribution of events: antiproton multiplicity vs. 
transverse radius for 10,000 simulated events for the scenario 
II. 



In this work, we have assumed that local equilibrium [1] 

holds, so that all variances and covariances are deter- 
mined by the equation of state and the conservation laws. 
When equilibrium is not in force, there are very few con- 
straints on correlations. We study the effect of incom- 
plete chemical equilibrium on the variances and covari- 
ances in |Kj]. While such results can be more realistic, 
they are also more model dependent. It will eventually 1 J 

be necessary to use models like RQMD, UrQMD, and 
VNI to study correlated observables. Such models per- 
haps offer the most realistic description of single par- 
ticle observables, but may require modification to ad- 
dress event-by-event correlations. For example, it has 
not been demonstrated that fluctuations in these models 
respect the fluctuation-dissipation theorem when one im- 
poses the appropriate boundary conditions and approxi- 
mate local equilbrium [ ]14j . Consequently, it is not clear 
that the dynamical correlations these models produce are 
realistic. 

On the experimental side, it is not clear how to extract 
event- by-event HBT radii in practice [Q . If experimen- 
tal fluctuations in HBT radii are too large, it will be 
necessary to develop a super-event analysis scheme f2l|| . 

We have shown that under the right circumstances a 
correlated analysis of observables such as antiproton yield 
and the pion Rt will augment the discovery potential for 
novel event classes. In particular, we stress the impor- 
tance of looking for systematic changes in correlations as 
functions of the centrality. A first study of fluctuations \J\ 

of the pion transverse momentum and the K/tt ratio by 
the NA49 experiment at the SPS focused on central col- 
lisions [J17J- Their study has similar motivations, such as 
testing the degree of equilibration (lqjn| and looking for 
critical fluctuations ||. Centrality studies would further ^ ' 

allow them to turn equilibrium "off" or "on." 
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